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Abstract 

Using the determinantal formula of Biane, Bougerol, and O'Connell, we give multitime joint 
probability densities to the noncolliding Brownian motion with drift, where the number of particles 
is finite. We study a special case such that the initial positions of particles are equidistant with 
a period a and the values of drift coefficients are well-ordered with a scale a. We show that, at 
each time t > 0, the single-time probability density of particle system is exactly transformed to the 
biorthogonal Stieltjes-Wigert matrix model in the Chern-Simons theory introduced by Dolivet and 
Tierz. Here one-parameter extensions (^-extensions) of the Stieltjes-Wigert polynomials, which are 
themselves g-extensions of the Hermite polynomials, play an essential role. The two parameters a 
and a of the process combined with time t are mapped to the parameters q and 8 of the biorthogonal 
polynomials. By the transformation of normalization factor of our probability density, the partition 
function of the Chern-Simons matrix model is readily obtained. We study the determinantal 
structure of the matrix model and prove that, at each time t > 0, the present noncolliding Brownian 
motion with drift is a determinantal point process, in the sense that any correlation function is 
given by a determinant governed by a single integral kernel called the correlation kernel. Using 
the obtained correlation kernel, we study time evolution of the noncolliding Brownian motion with 
drift. 
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I. INTRODUCTION 



Vicious walker models on lattices [1] and their continuum versions [2], which will be 
called noncolliding diffusion processes [3, 4], have been extensively studied in connection 
with the random matrix theory [5-14], the Tracy- Widom distributions and extreme value 
statistics [15-23], the enumerative combinatorics and representation theory [24-28], the 
Riemann-Hilbert problem [29-31], the renormalization group theory [32-34], the growth 
models [35-37], the quantum integrable systems [38-40], and others. Construction of non- 
colliding diffusion processes has been based on the determinantal formula for nonintersecting 
paths of Karlin-McGregor (KM) [41] and Lindstrom-Gessel-Viennot (LGV) [42, 43]. It can 
be regarded as a stochastic version of Slater's determinantal wave function for free fermions 
in quantum mechanics and it means that indistinguishability of particles (i.e. invariance of 
statistics under any exchange of paths at intersecting points on the spatio-temporal plane) 
should be assumed. In general, if each walker (diffusion particle) has different drift, then 
they are distinguishable from each other and the KM-LGV determinantal formula is not 
applicable. 

If the values of drift coefficients of many particles are well-ordered, however, Biane, 
Bougerol and O'Connell [44] proved that the following determinantal formula is valid. Let 
N = 2,3,..., and we consider an iV-particle system of Brownian motions in the one- 
dimensional real space R such that the j-th Brownian particle starts at Xj G R at time 
t — and it has a constant drift Vj per time, 1 < j ' < N. In other words, we consider an 
iV-dimensional vector-valued diffusion process 

B x ' u {t) = B x (t) + ut, t>0, (1) 

whose j-th component describes the j-th Brownian motion with a drift coefficient Uj, 

= Xj + Bj(t) + Vjt, t>0, l<j<N, (2) 

where Bj(t),t > 0, 1 < j < N, are independent one-dimensional standard Brownian motions 
all started at the origin. Let 

Wat = {x = (xi, x 2 ,..., x N ) 6l JV :i 1 <i 2 <---< x N }, (3) 
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which is called the Weyl chamber of type Ajv_i in the representation theory [45]. Biane, 
Bougerol and O'Connell (BBO) showed that [44], if the values of the components of initial 
configuration x — (x±, X2, ■ ■ ■ , %) and the drift vector v = (vi, v 2 , . . . , v N ) are in the same 
order, that is, 

x e Wjv, v e W N , (4) 

then the transition probability density of the drifted Brownian motions (1) conditioned never 
to collide with each other is given by 

det [e v ' Vk ] 

p»(t,y\x) = ^ m2/2 lJ ^T^ ^(t,y\x), yeW N , t>0, (5) 

l<j,k<N [ 

N 



where \v\ 2 = J2j=i v l anc ^ 



q N {t,y\x)= det p(t,yj\x k ) (6) 

l<3,k<N L 



with the transition probability density of B(t), 

p(t,y\x) = -l= e -^ 2 ^, t>0. (7) 

(We give a brief review of the BBO argument [44] in Appendix A.) 

In the formula (5), we note that even when some of Vj's in v = (yi, . . . , v^) coincide, the 
ratio of determinants deti< ? - j fc<jv[e I/j ' J ' fc ]/ deti<j tk < N [e L ' : > Xk ] can be interpreted using l'Hopital's 
rule. In particular, if we take the limit Vj — >■ 0, 1 < j < N, (5) is reduced to be 

p N (t,y\ x ) = f ^lq N (t,y\x), yeW N , t>0 (8) 
h N (x) 

with the Vandermonde determinant 

h N (x)= det [x k ,- 1 ]= TT (x k - Xj ). (9) 

Kj,k<N 3 J-- 1 - J 

It is considered as Doob's harmonic transform (/i-transform) [4] of the KM-LGV determinant 
gjv(t, y\x) by hjy, in which is harmonic in the sense that A?in(x) = ^2f = i d 2 hj^(x) /dx? = 
0. The above observation implies that the BBO formula (5) is a generalization of the h- 
transform of KM-LGV determinant (8). 

In (8), we can take the limit Xj — >■ 0, 1 < j < JV, which is denoted by x — > 0, and we 
obtain the probability density of particle positions y = (yi,y2, ■ ■ ■ ,Vn) of the noncolliding 
Brownian motions all started at the origin, 

PN (t,y\0) = c N (t)(h N (y)) 2 e-\y\ 2 l 2 \ y G W N , t>0 (10) 



with c N (t) = t-"* /{(2ir) N / 2 Uf =i r(J)}. The important fact is that (10) can be identified 
with the probability density of eigenvalues of random matrices (ordered in Wat) in the 
Gaussian unitary ensemble (GUE) with variance a 2 = t [46, 47]. When x ^ 0, (8) gives 
the probability density of eigenvalues {yj}jLi in the GOE-GUE two-matrix model studied 
in the high-energy physics [46], in which the hermitian random matrix is coupled with a 
real-symmetric random matrix with eigenvalues {xj}^ =1 [11, 48]. In this sense, the BBO 
formula (5) is expected to be related with some matrix-models which are generalizations of 
two-matrix models. (See [44, 49, 50] for the relations of (5) with the Duistermaat-Heckman 
measure and the Littelmann path model in representation theory.) 

On the other hand, the BBO formula (5) is regarded as a simplified version of the tran- 
sition probability density of the O'Connell process, which is a stochastic version of the 
quantum Toda lattice [39]. The reduction from the O'Connell process to the noncolliding 
Brownian motion with drift is called a combinatorial limit (or tropicalization) and the in- 
verse procedure is called a geometric lifting [49-52]. (See Appendix B.) In an earlier paper 
[51], we took the limit x — > with keeping v e W^r and derived a reciprocal time relation 
between the noncolliding Brownian motion with drift v started at and that without drift 
started at a configuration given by v. We would like to say that the BBO formula (5) is 
located in the high level in the random matrix theory and it also gives an introduction to 
mathematical models in higher levels [50, 52, 53]. 

In the present paper, we consider a special case such that the drift coefficients and initial 
positions are given by the following. Let n e N = {1, 2, ...,}, N — 2n — 1 and define 

/ N - 1 N- 1 N- 1 N-l\ 

P= (-— .-— + 1,.., -1,0,1. -—-1.—) 

= (-n + 1, -n + 2,--- ,n-2,ra- 1). (11) 

Then we put 

Xj = aPj = a(j -n) = a(j- ^p^J , 1 < 3 < N , ( 12 ) 

v i = °Pj = a ti -n) = a^j- —Y^j » l<j<N, (13) 

with positive constants a > 0, a > 0. Although this setting seems to be very special, 
we find that the obtained process can be transformed to a matrix model in Chern-Simons 
theory recently introduced by Dolivet and Tierz [54]. Precisely speaking, the normalization 
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factor of the single-time probability density in our process gives the inverse of partition 
function Z of their matrix model. Since here we study an interacting particle system (the 
noncolliding Brownian motion with drift) we can discuss not only the partition function 
but also correlation functions. Dolivet and Tierz shows that the ensemble of eigenvalues 
of their matrix model realizes the biorthogonal ensemble [55, 56] associated with the one- 
parameter extension [54] of the Stieltjes-Wigert polynomials [57, 58]. As an extension of 
their result, we will show in this paper that this ensemble is a determinantal (fermion) point 
process in the sense of [59, 60]. This implies that, at each time t > 0, the noncolliding 
Brownian motion started at the equidistant points (12) with the special drift (13) is also a 
determinantal point process, whose correlation kernel is expressed by using the biorthogonal 
Stieltjes-Wigert polynomials. We will report the time-evolution of the correlation functions 
of this noncolliding Brownian motion with drift. 

The paper is organized as follows. The probability densities of the noncolliding Brownian 
motion with drift and its transformations are given in Sec. II. The biorthogonal Stieltjes- 
Wigert polynomials are introduced in Sec. Ill and determinantal structure of the ensembles 
are studied there. In Sec. IV the main result is stated and analytic and numerical study of 
time-evolution of the present noncolliding Brownian motion with drift using the correlation 
kernel is reported. Concluding remarks will be given in Sec.V. Appendices are prepared for 
the BBO argument, the O'Connell process, and the geometric Brownian motions, which are 
related with the present process. 



II. PROBABILITY DENSITIES OF THE PARTICLE SYSTEMS 

A. Multitime Joint Probability Density of the Noncolliding Brownian Motion with 
Drift 

In a previous paper [51], we considered an A^-particle system of noncolliding Brownian 
motion started at x — (x 1 ,x 2 , ■ ■ ■ ,xn) G Wjv = {x G M N : x\ < x 2 < ■ ■ ■ < x N } with a 
drift vector v = (z/i, z/ 2 , . . . , v N ) satisfying u G Wat. For any M G N and an arbitrary M 
sequence of times < t\ < t 2 < ■ ■ ■ < tu < oo, the multitime joint probability density of 
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the process is given by 

P %(t 1 ,xW;...',t M ,xM\ x ) 



e -t M \v?li det 

l<j,k<N 



(M) 
e V 3 X k 



M-l 



m=l 



det f. 

where particle configurations at each time t m is denoted by x^ = (x^, . . . ,xffl) G 
W^, 1 < m < M. Here is the KM-LGV determinant given by (6). From this gen- 
eral formula, we can obtain the following multitime joint probability density for the special 
case (12) and (13). 



Proposition 1 Consider the noncolliding Brownian motion with iV particles started at 
the equidistant points (12) at time t — with the drift having the coefficients (13). For 
an arbitrary M G N and an arbitrary sequence of times < t\ < t 2 < ■ • • < t M < oo, the 
multitime joint probability density of the process is given by 

N 

p N (h, • • • ; t M , * (M) ) = c N (a, a, t u t M ) JJ( e -r )-(^-D/2 



i=i 



x Y[ (e^ M) -e^' M) ) Y[q N (t m+1 -t m ,x^\x^) 

l<j<k<N m=l 
N N 



M-l 



1 



with 

c N (a,a,ti,t M ) = 
It is also written as 

p N (ti,x^;...;t M ,x^) 



Yl^ e naa _ l)N -n 



l<j<k<N 



exp { -^N(N 2 - 1) (a 2 t M - 2aa + ^- 



2 sinh 



= c N (a,a,ti,t M ) I| 

l<j<k<N 
N 

xjipttux^o) n 

j=l l<j<k<N 



M-l 



(15) 



(16) 



n^(wi-^^ (m+i v (m) ) 



m=l 



^ . , a \ X k ~ X i ) 

2 sinh — 3 — 

2ti 



(17) 



Proof If we set the drift vector as (13) with a positive constant a > 0, we obtain 

AT 

-(7V-l)/2 -Q 



det e"'' Xfc = det 

l<j,k<N l<j,k<N 



( e <rx k y-l-{N-l)/2 | = Jj^ e <rz^-(iV-l)/2 jj (g<rs fc _ e <rx^ _ 
£=1 l<j<k<N 
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We find that the above is equal to 



n 

l<j<k<N 



2 sinh 



a(x k - xj 



Note that 



A? 



\u\ =a\p\ = a 2^ I / 



Then the multitime joint probability density of the process with drift coefficients (13) is 
given by 



Pn 



°P(t 1 ,zV;...;t M ,tM\x) 



N 



= e 



n 



iV(iV a -l)<T 3 t M /24TT/.<TxJ M) \-(JV-l)/2 



n 



(M) (M) 



l<j<k<N 



M-l 



x n^(wi-^^ (m+i) k w )- 



^(tijaj^lx) 



m=l 



= e 



-Af(Af 2 ^l)cr 2 t M /24 



n 



J-J( era ,)-(K-l)/2 J-j ( e ™»_ e ™,) 



2 sinh 



l<j<k<N 



M-l 



9jv(*i,* (1) I 



m=l 



n 

l<j<k<N L 



2 sinh 



(18) 



Then we assume that the initial configuration is given by the equidistant points (12) with 
a > 0. We can see that in this setting 



q N (t 1} x {1) \ap) 



D -N(N 2 -l)a 2 /24ti e 



\X^\ 2 /2t 1 



AT AT 

= e -A^ 2 -l)« 2 /24 tl JJp^^WjO) JJ( e ax^/ti)-(^-l)/2 J] (gOS^/fi _ e «f >/*i; 

l<j<k<N 
I (1) (1) 

a(4 - x} 



j'=i 

Af 



<=1 



-Af(Af 2 -l)a 2 /24t 



l<j<k<N 



2 sinh 



2ti 



Moreover, we can show 



v 



Y^^a(e-(N+l)/2)y(N-l)/2 = ^ 



1=1 
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and 



JV-1 



JJ | e ^(fc-(iV+l)/2) _ e aa(j-(AT+l)/2)| = e -JV(JV 2 -l)oa/12 JJ ^naa _ ^AT-n 
l<j<k<N n=l 

Therefore, we obtain the expressions (15) and (17) with (16) for p N (- ■■)= p a N p (- ■ ■ \ap). 



B. Transformation to the System Associated with Geometric Brownian Motions 



Let M + = (0, oo). We change the variables as x" e R — > y E R+ by 



e ax > 



(m) 



,M 



x 



(m) 



- lW m) , 1 < m < M, l<j<N, 



(19) 



and put 



M 



M 



p N (h, ■■■;t M , ajW) H <te (m) = p*(fi, • • • ; *m, ?/ M) ) H d f 



(m) 



m=l 



m=l 



with rfa;( m ) = nf=i ^4 m) , ^ (m) = Fill d Vj m \ l<m<M. Then we obtain 



N 



PN(h, y {1) ; ■■■;t M , y {M) ) = c N (a, a, t u t M ) \\{y 



(Af)x_(AT-l)/2 



£=1 



M-l 



J'=l 



with 



and 



where 



P geo (t,y\x) 



yaylut 



exp 



l<j<k<N 



(My MY 

2aH 



q g ™(t,y\x) = det p geo (£, 2/j-|ar fc ) 

l<3,k<N L 



(20) 



x n (yr-ynll^Wi-^y^li/^) 

i<j<k<N "i=i 

N N 

xn^^^fii)!!^)"^- 1 ^ 072 n - (^ (i v (ti) ] (2i) 



(22) 



(23) 



(24) 



If we see (22) as a function of y, it is considered as the probability density of the log-normal 
distribution with parameters a\/i and logrr. As explained in Appendix C, it gives the 
transition probability density from x to y in time duration t > of the geometric Brownian 
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motion X(t),t > 0, which is given by an exponential of one-dimensional standard Brownian 
motion B(t); 

X(t) =xexp(aB(t)), t > (25) 

with the initial value x = X(0) and with the parameter a > called the percentage volatility 
[61]. 

It should be noted that, though (23) is the KM-LGV determinant of p sco, s, (21) is different 
from the transition probability density of 'noncolliding geometric Brownian motion'. (See 
Appendix D for detail.) 



C. Single-Time Probability Density and Transformation to Biorthogonal Ensem- 



ble 



From now on we consider only a single-time distribution. By setting M = 1 and simpli- 
fying the notations as t — t\, y — y^\ (21) gives 



PN(t,y) = c N (a,a,t,t)Y\ 



j'=i 



1 



y. 



-(AT-l)0(t)/2-(JV+l)/2 



exp 



2aH 



x n iy k -vM {t) -y?\ 

l<j<k<N 



(26) 



If we change the variables as yj — > Zj by 



exp 



Vj = zj, l<j<N, 



(27) 



we obtain the probability density of the distribution z = (z 1: z 2: ■ ■ . , z N ) G = {z e M. N : 
< Z\ < z 2 < ■ ■ ■ < zjy} (the Weyl chamber of type Cat), which is given by 

N 



P N (t,z) = C N (a,a,t)l[w(z J ) J] {(z k - Zj )(z^ - zf } )}, 



(28) 



l<j<k<N 



where 



C N (a,a,t) 



lln=l 



x exp 



N 



12 



--UN + 1){2N + l)a 2 t + 2{N 2 - l)aa + (N — 1)(2N - 1) — 



and 



W(z) = l- e -^ogzf with p = _J_ 



(29) 



(30) 
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The probability density of the GUE (10) is proportional to a square of the Vandermonde 
determinant, while (28) is proportional to the product of 1in(z) and h,N({Zj}). The ensemble 
with the probability density in this form (28) is called the biorthogonal ensemble and studied 
in [55, 56]. The special case with the weight function (30) was studied in the name of the 
biorthogonal Stieltjes- Wigert matrix model for the Chern-Simons theory by Dolivet and Tierz 
[54, 62], since (30) is the weight function for the Stieltjes- Wigert orthogonal polynomials 
[57]. In particular, when t — to = a/a, 6{t) becomes unity and the system is reduced to the 
Stieltjes- Wigert matrix model studied by Tierz [63]. See Remark 3 in Sec. Ill C. (Since the 
partition function Z is given by a hermitian-matrix integral as Eq.(1.6) in [54], it is called 
the matrix model. The partition function Z is also written as the integral of weights of real 
eigenvalues as Eq.(31) below, then the matrix model is identified with a statistical ensemble 
of points on R.) 

Remark 1 In the theory of Chern-Simons matrix models [54, 62-65], the main quantity 
to be calculated is the partition function, which will be given by 

r N 

Z= Hw( Zj ) J] {(z k -z 3 )(z e k (t) -z 6 ^)\dz. (31) 

Jw S j=l l<j<k<N 1 J 

In the present setting, through the relations (24) and (30), Z is a function of a, a and t. As a 
matter of course, in order to identify (31) with the partition function studied in [54], we have 
to rewrite these parameters by using the proper parameters in the Chern-Simons theory, but 
the partition function is essentially given by (31). Since P/v(t, z) given by (28) is a probability 
density, it is normalized, and then (31) is equal to the inverse of Cjy given by (29). Here we 
would like to put emphasis on the fact that, since in the present work the "matrix model" 
is realized as a transform of the stochastic process (the noncolliding Brownian motion with 
drift), Cn has been automatically obtained by the transformation from the normalization 
factor c N given by (16), and cjy is readily obtained by just putting the proper conditions 
on the initial configuration (12) and the drift vector (13) as demonstrated in the proof 
of Proposition 1. In [54] the partition function was calculated by using the biorthogonal 
Stieltjes- Wigert polynomials. We do not need them to evaluate Cn = 1/Z, but we will also 
introduce them in the next section in order to discuss the correlation functions of the matrix 
model and our stochastic process. 
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III. TIME-DEPENDENT BIORTHOGONAL STIELTJES-WIGERT ENSEMBLE 
AS A FAMILY OF DETERMINANTAL POINT PROCESSES 



A. Biorthonormal Stieltjes-Wigert Polynomials 

As functions of q G (0, 1), we consider the g-extension of the Pochhammer symbol as 
(a;q) = 1, 

(a; q) n = (1 - a)(l - aq)(l - aq 2 ) ... (1 - aq 71 ' 1 ), n eN, 



for aeR, and the g-binomial coefficients defined by 

(?;?)n n 



n 

I 



l,2,...,n-l, 



and 



n 




n 





i 


n 



= 1. 



Let 



: CXp 



(In z) 



(32) 



y/2ir\\nq\ ~" r I 2 1 In g | 
The following two series of functions were introduced in Appendix A. 2 in [54], for 9 G 
(0,oo),ge (0,1), 



, Hn- n\ n( ne + 1 / 2 )/ 2 n 

r.(« ; «, 9 ) = (-i)»V(!^» £(-!)< 

{q ,q )n e=0 



n 
I 



(nfl+l/2)/2 » 

iZ B (x;g,g) = (-ir g /r -^ ^(-l) < 

\/{q\q)n e=0 



n 

I 



q e(e(e+i)+(i-d)+i)/2 x e^ 



(33) 



n G M = {0, 1, 2, . . . }, where T„(x; 6 1 , g) is a polynomial of a; 6 of order n and i?„(x; 9, q) is a 
polynomial of x of order n, n G N . 

Proposition 2 The polynomials {T n (x; 9, q), R n (x; 9, q)} n eN satisfy the following orthonor- 
mality relation with respect to the weight function (32), 



POD 

/ T n (x;9,q)R m (x;9,q)w(x;q)dx = 5 nm , n,mGN . 
Jo 



(34) 
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Remark 2 We can see that 

UmT n (x;9,q) = \im R n (x; 9 , q) = p n (x;q), n G N , (35) 
where p n (x; g),n G N are the Stieltjes-Wigert polynomials given by [57] 

q i2 (-q^xY, nG N . (36) 



(2n+l)/4 n 

p B (x;9) = (-ir^ 7 ==E 
V(9;?)» £=o 



(Note that the definition (36) is slightly different from that given in [58].) Dolivet and 
Tierz [54] derived {T n (x; 6, q), R n {x\ 9, q)} n m by taking a limit a — > oo of the g-Konhauser 
polynomials, which has a parameter a in addition to 9 and q. The weight function, with 
which the g-Konhauser polynomials make orthogonality relations, is called the g-Laguerre 
measure [66]. We note that the a — > oo limit of the g-Laguerre measure is different from the 
present weight function w(z; q) given by (32). (It is related with the weight function used in 
[58] to define the Stieltjes-Wigert polynomials.) We give a proof of the orthonormality (34) 
with respect to (32) below. It implies that as well as the Stieltjes-Wigert polynomials (36), 
the biorthonormal polynomials {T n (x; 9, q), R n (x; 9, g)} n eN are indeterminate and there are 
many different weight functions (the Stieltjes moment problem, see [63, 65]). 

Proof of Proposition 2 It is enough to prove the following, 

poo 

I n ,m = / x m T n (x; 9, q)w(x : q)dx = t n {9, q)5 nm , < m < n, (37) 
Jo 

poo 

J n ,m = / x m R n (x; 9, q)w(x : q)dx = r n (9, q)5 nrn , < m < n, (38) 
Jo 

where 



tn(9,q) = Vi^Uq'^ M'-W 

rn (0, q ) = (gVk g-M+i/^^^-i)^ (39) 



It is easy to confirm 



poo 

/ x n w(x;q)dx = q- (n+1)2/2 , n G N (40) 
Jo 

for the weight function (32). Then by definition of the polynomials (33), 

J — A V lg i q >A n {n-m)te -(m+l) 2 /2 /^x 

V- U "Z^ (a d -a d )p ' ^ ' 

e=o ^ ' y > l 

n ( n -n6. n 9\ 

Jn,m-d n ^ {qe ^ e)( q q , (42) 

12 



with 

J(a-a) ( ne + 1 /2)/2 (nS+l/2)/2 
dn = (-1)" V(g '?y „ , rfn = (-l)^ 



The g-derivative of order n of a function /(x) is defined as [66] 



(43) 



* M ^ (l_ g )n x nZ. M« ^ 



Then 



i n („-n9.„9\ 

Let f(x) = x m qr-( m + 1 ) 2 / 2 in (45) and compare the result with (41). We obtain 

I n ,m = 4(1 - q -<>y q -(™+V 2 /* D n q _ e x m \ x=i . (46) 

It implies that J„ im = for < m < n. For n = m, by direct calculation, we see I n ^ n = 
t n (9,q). Then (37) is proved. For the proof of (38), put f(x) = x ™ q -^+if/2 in (45). We 
obtain 

J n ,m = 4(1 " q-y q ^+V 2 l* D n q - e x m \ x=1 (47) 
and (38) is derived in the same way as (37). Then the proof is completed. 1 

B. Time-Dependent Stieltjes-Wigert Ensemble 

For a, a > 0, let 

= 6{t) = —, q = q(t) = e~ u2 \ t > 0. (48) 
at 

Then, by the equality between the Vandermonde determinant and the product of differences 
(9) and by the multi-linearity of determinant, the probability density (28) of z is written as 
follows, 

N 

P*(t, z) = Hw( Zj ; q(t)) det b-ife 6(t), q(t))] det \R m -i(zf, 9(t), q(t))] . (49) 

i<J,k<N L J l<t,m<N L J 

i=i 

It is a one-parameter family with parameter t > of the biorthogonal Stieltjes-Wigert 
ensembles. We call it the time- dependent biorthogonal Stieltjes-Wigert ensemble. 
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C. One- Parameter Family of Determinantal Point Processes Parameterized by 
Time 



Let C (R+) be the set of all continuous real- valued functions with compact support on 
R+. For / e C (R+) and ft £ 1, at each fixed time t > 0, the generating function of 
correlation functions is given by the following Laplace transform of PN{t, z), 



g N (t, X ) = [ e«^™P N (t,z)dz 

i r N 

= m N U.0--x(zj))PN(t,z)dz, 



(50) 



where we put 



X (z) = 1 - e*'«. (51) 
If we write Zn> = (zi, Z2, ■ ■ ■ , zw), N' < N, the binomial expansion of the integrand of (50) 



gives 



N N' 

g N (t, X ) = l+J2(-l) N '^ I Ni PN(t,z N ,)H{x(z k )dz k }, (52) 

N'=l ' J R +' k=l 



where p^it^ZNi) is the iV'-point correlation function at time t > 0, 



N 



p N (t, z N> ) = W ± WY _ J^ N _ Ni P N (t, z) IJ dz jt 1 < N> < N. 



(53) 



By the determinantal expression (49) of the single-time probability density -P/v(t, z), (50) 
is written as 



GN(t,x) 



det 

l<j,k<N 



T^zf, 9(t), q(t))w(zj; q(t))(l - X {zj))] 



x det 

K£,m<N 



R m ^(zf,e(t),q(t)) 



dz, 



(54) 



where multi-linearity of determinants is used. For square integrable functions gj, gj,l < j < 
N on A C R, the identity 



. . . det 

N\ J a n i<j,k<N 



9j{xk) 



det 

J K£,m<N 



ge{x m ) 



dx = det 

l<j,k<N 



gj(x)g k (x)dx 



(55) 



is proved, which is called the Andreief identity. Then, if we set 



Bj(z) = T^z; 6(t), q(t))y/w(z; q(t)) X (z), 
Cj(z) = Rj.^z; 9(t), q(t))y/w(z; q(t)), j G N , 



(56) 
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we have the determinantal expression 



GN(t,x) = det 

l<j,k<N 



poo 

5 jk - / Bj(z)C k (z)dz 
Jo 



(57) 



where we have used the orthonormality (34) proved in Proposition 2. 

Then, by Fredholm's expansion-formula of determinant and by cyclic property and multi- 
linearity of determinants, we have the equalities 



N 

N'=l JK + ~ 



N 



e=i 



dzMi 



N 



= i + E<-i: 



N' 



, 1 



JV'=1 



N'\ i<i 



N' 

d &\KN(t;z j7 z k )} TT \x(z k )dz k \ (58) 

j,k<N' ) 
k=l 



with the integral kernel for (x, y) G M? + with t > 



JV-l 



K N (t;x,y) = y/w{x; q(t))w(y; q(t)) ^ ^(x; 9(t), q(t))Tj(y; 0(t), q(t)). 

3=0 

The rhs of (58) is the definition of the Fredholm determinant, which is denoted by 



GnQ, x) = Det a . )W6M+ 5{x-y)- K N (t; x, y)x{y) 



(59) 



(60) 



Comparing (58) with (52), we can conclude that, at each fixed time t > 0, for any 1 < N' < 
N, the A^'-point correlation function is given by the determinant in the form 



p N (t,z N >)= det K N (t;Zj,z k ) , 1 < N' < N, t > 0. 
In particular, the particle density is given by the one-point correlation function as 



(61) 



p N (t;z) = K N (t,z,z), z > 0, £ > 0. 



(62) 



The ensemble of points such that any correlation function is given by a determinant (and 
thus the generating function of correlation functions is given by a Fredholm determinant) is 
called the determinantal (or fermion) point process [59, 60]. The integral kernel is called 
the correlation kernel. 

Remark 3 For each a > 0, a > 0, there is a special time 



(63) 
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at which 



a 2 

0(*o) = 1, v% = aa = —, q(to) = e~ aa = q . (64) 

to 

At t = t , the probability density (49) becomes 

N 

P N (t -z)=c N (q )l[w(z :j] q ) J] (z k - z 3 f = F$(z) (65) 

j=l l<j<k<N 

with 

N(W 2 -l)/6 

c N {qo) = ^n-!, -• ( 66 ) 

llj=i («b; 9b)j 

Then the system is reduced to the Stieltjes-Wigert ensemble studied by [63], which is also a 
determinantal point process with the correlation kernel 

AT-1 

K n(x,v) = \/w(x; q )w(y; q ) p 3 (x; q )Pj(y, q ), (x,y) G R+, (67) 

3=0 

as derived from (59) by (35). For the Stieltjes-Wigert polynomial p n given by (36), the 
three-term recurrence relation is given by 



Pn{x; q) = Pn-i{x; q)-q — , Pn-2{x; q), n = 2, 3, . . . . (68) 

a/1 - q n V 1 - g n 

We can obtain the Christoffel-Darboux formula for them, 

^ r v , x _ y/l-q N p N (x;q)p N - 1 (y; q) - Pn-i(x; q)p N (y; q) 
2_^Pn\x,q)p n {y,q) - , 

n=0 9 X ^ 

for x ^ y, 

x- 1 R ~ q N f . 

^p„(:r) 2 = — — ^p' N (x;q)p N _ 1 (x;q) - p' N _ 1 (x;q)p N (x;q)j. (69) 

n=0 ^ 

Then the correlation kernel (67) is rewritten as 

7^0/ \ V 1 ~ Vo n w ^PN(x;q )pN-i(y,qo) -PN-i(x;q )pN(y,qo) 

K q N(x,y) = m — y/w[x\q )w{y\q Q ) , 

% x y 

for x ± y, (70) 

R _ g N , . 

k n(x,x) = Tr ^w(x;q )\p' N (x;q )pN-i(x;qo) - p' N ^(x; q )p N (x; q )\ ■ (71) 

Remark 4 Consider the system again at the special time t — to, but here we write qo = q 
for simplicity. Let 

3/2 

z^q-^yW^ + l- 1 ' 2 <=► z, = "J^ , l<j<N. (72) 
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Then we can show that 

N N 

lim F*(z) J] dzj = Pk(z) J] dzj (73) 

Q ^ 3=1 3=1 



with 



-N 2 /2 v 



P ^ ^^m W\MJ^ 56W " (74) 

By (64), q = q = e~ au — >■ 1 as acr — >■ 0. Note that (74) is equal to the special case with the 
variance a 2 = 1/2 of the probability density of eigenvalue distribution of the GUE, 

N 2 N 

P gue (a) = * TT^W TT ( A A )2 AGW Wl (75) 

(see (10) in Sec. I.) Define 
K q N (x,y) 

= K% (q-^x v^T^) + g- 1 / 2 , g" 3 ' V^O^) + <T 1/2 ) <T 3/2 v^O^) , (76) 
where K^(-, •) is given by (67) with q = q. The following asymptotics is established 



lim ^^ w/2 q' n/2 - 1/4 Pn (q~ 3/2 VW^7)x + q^ 2 ) = H n (x), n G N , (77) 



{(l"?)/2}' 

where H n (x),n G No are the Hermite polynomials, 



fc=0 

Then, as expected, we can prove the following convergence, 

N-l 

Um K q N (x,y) = V ^(x)^(y), (78) 

q— >1 ' 

j=0 

where y?j(:r) is the Hermite function, <£>j(x) — e~ x2 / 2 Hj (x) / a/ 2-? j! y^, j G No- 

Combining the results in Remarks 3 and 4, we can say that in the double limit 9—^1 and 
q — >■ 1, the biorthogonal Stieltjes-Wigert ensemble is reduced to the GUE ensemble, which 
is a determinantal point process with the Hermite kernel (78) under an appropriate scaling 
of variables (72). 
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IV. TIME EVOLUTION OF THE NONCOLLIDING BROWNIAN MOTION 
WITH DRIFT 



A. Main Result 

By the sequence of transformations (19) and (27), each variable Zj of the time-dependent 
Stieltjes-Wigert ensemble (49) is related with the variable Xj = xf\t = t x for the original 
noncolliding Brownian motion with drift given in Proposition 1 in Sec. II A. That is, 



Zj = exp 



exp 



ax j + -{(N-l)e(t) + (N + l)}aH 



axj + i(JV - l)aa + ]-{N + l)a 2 t 

Z Z 



(79) 



1 < j < N, where the relation (24) was used. Since dzj = azjdxj, 1 < j < A, we arrive at 
the following main result of the present paper, where we have used the cyclic property of 
determinant. 

Theorem 3 Consider the noncolliding Brownian motion with iV particles started at the 
equidistant points (12) at time t = with the drift coefficients (13). At each time t > 0, 
the particle configuration {Xj(t)}^ =1 is the determinantal point process in the sense that 
any A'-point correlation function of x N / = (x x , x 2 , ■ ■ ■ , Xw), 1 < A' < A, is given by 
determinant 



PN(t, Xn') — det 

l<j,k<N' 



K N (t;Xj,x k ) 



(80) 



Here the correlation kernel K N is given by 

K N (t;x,y) = K N {t-e^ N -^' 2 ^ N+1 ^ 2t ' 2 ^y+(N-i)aa/2 H N + i) 



2v 



xae *(x+y)/2+(N-l)aa/2+(N+l)* 2 t/2^ £ 



(81) 



where K N (t;x,y) is given by (59) with (24). 



B. Time-Evolution of Particle Density 

The particle density function for the A-particle process is given by the one-point (A' = 1) 
correlation function, which is a 'diagonal value' of the correlation kernel, 



p N (t,x) = K N (t;x,x), iGK, t > 0. 



(82) 
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FIG. 1: The particle-density profile /5jv(i, x) of the noncolliding Brownian motion with drift at 
time t = 1 for N = 15, where a = a = 1. The number of nodes is equal to N = 15. The oscillatory 
behavior represents the lattice structure in one dimension caused by the present setting of the 
equidistant initial configuration and well-ordered drift coefficients. 

Figure 1 shows (82) for iV = 15 and t — 1, a — a — 1, i.e., 9 — 1 and q = e~ l ~ 0.37. The 
oscillatory behavior of density profile is observed, which was already reported in [65]. In 
the present work, we put the equidistant initial configuration (12) and the drift coefficients 
which are regularly ordered as (13), and then the system has a lattice structure in one 
dimension represented by this oscillatory behavior. 

We compare the dependence on time t and particle-number N of density functions be- 
tween the present process, pAr(t,x), and the noncolliding Brownian motion without drift 
started at 0, denoted by P7v(t, x). As reviewed in Sec. 5.1 in [3], the latter is given by 

n\ Vn ( -|= ] 1 - y/N(N + [ — ] r V-1 



X 



(83) 



x G R, t > 0, where tpj is the Hermite function, and it behaves asymptotically in iV — > oo as 
p N (t, x) i 



-J^Ntf-x^ 



(84) 



if \x\ < 2VNt, 
0, if |x| > 2VWt. 

Figure 2 shows (83) at time t = 1 for N = 1,5,9,13,17 (from inner to outer). We can 
observe oscillatory behavior also in this process (the Dyson's Brownian motion model with 
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FIG. 2: The particle density functions pisr(t, x) are plotted at time t = 1 for N = 1, 5, 9, 13, 17, from 
inner to outer, for the noncolliding Brownian motion without drift started at 0. As N becomes 
large, the profile becomes a semicircle with radius oc y/N. 
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FIG. 3: The particle density functions pjsr(t,x) are plotted at time t = 1 for N = 1,5,9,13,17, 
from inner to outer, for the present noncolliding Brownian motion with drift (a = a = 1). As N 
becomes large, the support of the profile becomes wider proportionally to N. Then the height of 
the profile is constant and the lattice structure with N nodes is maintained. 

(3 = 2 started at 0), but the support of density function of iV particles can expand only 
in the order a/ Nt. Therefore, as shown in Fig. 2, the 'wave length of oscillation' becomes 
smaller as \ft/N — > in — » oo for each time < t < oo and we will have Wigner's 
semicircle law (84) for the particle-density profile asymptotically in A^ — > oo. 
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FIG. 4: The iV-dependence of width of density profile A at time t = 1.5 is shown for a = a = 1 
(dotted by black circles), a = a = 0.5 (by black squares), a = a = 0.25 (by black triangles), and 
for a = a = (by white circles). For large values of N, the dependence is well described by a 
line A = c\N + c 2 for the three cases with a = a > 0; (ci,c 2 ) = (2.50,6.29) for a = a = 1, 
(ci,c 2 ) = (1.25,9.70) for a = a = 0.5, (ci,c 2 ) = (0.627,15.9) for a = a = 0.25, respectively. The 
white circles for the noncolliding Brownian motion without drift started at approaches to the 
result of Wigner's semicircle law for the GUE described by a curve A = 4V1.5JV ~ A.9VN. 

On the other hand, Fig.3 shows the particle density functions (82) of the present noncol- 
liding Brownian motion with drift with t = a = a = 1 for N = 1,5, 9, 13, 17 (from inner to 
outer) . 

As pointed out by de Haro and Tierz for the case with 9 — l,q < 1 [65], the width 
of support of profile increases proportionally to N, instead of \J~N . In order to see N- 
dependence numerically, here we define the width A(t) of density profile at a given time 
t > as the length of interval of x in which the value of p)v ( or Pn) is greater than 
e = 0.001. Figure 4 shows the iV-dependence of the width at time t = 1.5, A = A(1.5), 
for the cases a = a = 1,0.5,0.25 and for the case a = a = (the noncolliding Brownian 
motion without drift started at 0). By (48), for the first three cases, 6 = 1.5 -1 ~ 0.67 in 




21 



common and q ~ 0.22,0.69, and 0.91, respectively. In the cases with a = <r>0(g<l 
and 9 ^ 1 in general), A increases linearly in N, while it seems to approach to the curve 
A = A\/l.5N ~ 4.9v / iV in the case a = a = as expected (Wigner's semicircle law). 

Since we put drifts as (13) to particles, t-dependence of the support-width is t (drifted), 
instead of y/i (diffusive). In summary, we conjecture that, at each time < t < oo, 

width of support of f>N(t, x) oc Nt as N — > oo. (85) 

Since N particles exist on the support with width oc N, the average height of density profile 
becomes independent of N (oc 1/t in time t) and the lattice structure (the oscillatory behav- 
ior with iV nodes) will not become indistinct even if N is large. Note that by noncolliding 
condition the particles will show the fermionic (exclusive) behavior. It is in contrast with 
Wigner's semicircle law (84), in which the height of profile increases as N becomes large; 
p N (t,0) oc y/WJl. 



V. CONCLUDING REMARKS 



In the present paper, we study the noncolliding Brownian motion started at N equidistant 
points (12) with a period a for which drift coefficients are chosen as (13) with a scale a of 
values. If we take the double limit a — > and a — > 0, the process is reduced to the 
noncolliding Brownian motion without drift started at 0. It is well-known that, in this limit, 
the particle-position distribution is equivalent with the eigenvalue distribution of random 
matrices in GUE with variance t and each time t > it is a determinantal point process 
with the correlation kernel expressed by the Hermite polynomials [46, 47]. We have shown 
that, for any choice of positive values of parameters (a, a), the determinantal structure of 
particle distribution is maintained, such that a pair of parameters (9, q) are determined at 
each time t > by 

q = q(t) = e-° H , (86) 
9 = 9(t) = ±, (87) 

and the correlation kernel Kn is expressed by the (9, g)-extensions of the Hermite polyno- 
mials, which are the biorthogonal Stieltjes-Wigert polynomials introduced by Dolivet and 
Tierz [54]. 
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As shown by (86), if the system does not have any drift term, a = 0, then q=l, while if 
a > 0, then q — 1 once at t — and g decreases monotonically in time t > 0. Introduction 
of drifts into the system is essential for the present g-extension and the derivation of value q 
from 1 measures the time t. In other words, the present stochastic process is a system with 
a time-developing g-parameter (86). By (87), we see that for each setting of (a, a), there is 
a unique critical time to = a/a at which 9 (to) = 1. 

Generalization of Theorem 3 for time correlation functions will be an interesting future 
problem. The Eynard-Mehta-type correlation kernel [3, 67-69] shall be discussed. 

Invalidity of Wigner's semicircle law for q < 1 is an interesting phenomenon, which was 
first observed by de Haro and Tierz in the case 9 = 1 [65] and is also reported for 9 ^ 1 
in Sec. IV B in the present paper. Asymptotic analysis of correlation kernel (81) in N — > oo 
will be studied so that the conjecture (85) is proved. 

It has been shown here that the present noncolliding Brownian motion with drift is exactly 
transformed to the biorthogonal Stieltjes-Wigert matrix model studied by Dolivet and Tierz 
[54]. We expect further connections in mathematics and physics between nonequilibrium 
statistical mechanics and the Chern-Simons theory in high-energy physics. 
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Appendix A: Drift Transform and the BBO Formula (5) 

First consider a one-dimensional standard Brownian motion without drift started at the 
origin, B(t),t > 0; 5(0) = 0. The probability density at time t > is given by u(t,x) = 
e -x 2 /2t i yj2nt. It solves the diffusion equation du(t,x)/dt = (l/2)d 2 u(t, x)/dx 2 and we can 
see that u(0,x) = \im t ^ u(t, x) = 5(x). Let i/Gl and consider the diffusion equation with 
a drift term 
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A solution is given by 



u u (t,x) = u(t,x + ist) 



1 expi- — (x + ut) 2 }. (A2) 



V2^t I 2t 

The mean of the Gaussian distribution is shifted by —id, that is, —v gives a drift velocity. 
We note that (A2) is written as u u (t,x) = e~ vx ~ v2t l 2 u{t,x). It is considered as follows (see, 
for instance, [70]); u u (t,x) is obtained from u(t, x) by the drift transform, 

, 2 
2 



u{t,x) — > exp I — vx — —t ) u{t,x). (A3) 



The transition probability density from x to y with time duration t is then obtained by a 
shift of spatial coordinate 

p u (t,y\x) = u u (t,x-y) 

= exp |z/(y - x) - yij p(f, y|x), (A4) 

where p(t,y\x) is given by (7). Note that (A4) satisfies the diffusion equation (Al) as a 
backward Kolmogorov equation. 

Let N — 2,3, ... . Consider the vicious Brownian motion defined as 

A^-particle system of Brownian motion killed when they collide 
= A^-dimensional Brownian motion in the Weyl chamber Wat with absorbing walls. 

The transition probability density from a; 6 Wjv to y 6 Waj with time duration t > is 
given by the KM-LGV determinant gjv(t, y\x), (6). It is a unique solution of the differential 



equation 



Q I N Q2 

(t,x) = -Au(t,x), A = ^2 o3> (A5) 



j=l 3 



satisfying the boundary condition u(t,x) = at x e dW N , and the initial condition 
u(0, x) = 5(x-y) = Ylj=i^j ~ Vi)- 

Now we consider the vicious Brownian motion problem with drift. For v = (vi, . . . , v N ) G 
M. N , we want to solve 

d 1 

—u(t,x) = -Au(t,x) + h>-Vu(t,x) (A6) 
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with the conditions 



u(t,x)>0, xeW N ,t>0, 
u(t, x) = at x e dW N , t > 0, 
u(0, x) = 5(x — y). 



The solution is given by the drift transform of the KM-LGV determinant (6), 
q%(t, y\x) = exp ji/ • (y - x) - ^-f J 9jv(*> y\x) 



exp < v ■ (y — x) 



l/ -t \ det 

2 | l<j,k<N 



P{t,yj\x k ) 



The important fact is that 



q%(t,y\x) ± det p"'(*» S/jl 3 *) 

l<j,k<N . 

That is, it is not the KM-LGV determinant of the drift transform of p's. Actually 
see 



A? 



g£(t,y|a:) = e~ ux ^ sgn(a) J] J]p(t, - 

ueSff k=i j=i 

N N 



o-eSjv fc=i i=i 

where <Sjv is the set of all permutations {<r} of iV items. 
The survival probability at time t > will be given by 

M%(t,x) = / q%(t,y\x)dy 



N . N 

aeSjv k=l Jwn j=l 



It was claimed in [44] that 



then 



if v e Wat and x G Wjv 



lim / IT tfi (t, yj \x a{j) )dy = 1, Vcr e 5 

That is, for any initial configuration a? e K^, if (All) is satisfied, 

lim B x ' u (t) e W N with probab. 1. 
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It is a matter of course, since we put the drift vector v G Wat, the position vector of 
Brownian motions should be in Wat if we wait for sufficiently long term, independently of 
initial configuration. Then we obtain [44] 

N 

]im M%(t,x) = e u x V sgn(a) TTe^w 

ct£5jv k=l 

= e vx det \e v ' Xk \. (A14) 

l<j,k<N L J 

Remark 5 The following argument is found in [50]. Let Sj^(x) be the survival probability 
for the noncolliding BM with drift v G Wat started at x G Wat. Then, it is a stationary 
solution of the diffusion equation 

= ^AS%(x) + v ■ VS%(x), (A15) 



satisfying the conditions 



S%(x) = at x G dW N , 

Hm Sj?(as) = 1, (A16) 



where 

a; — )• oo -<=>- — Xj — > oo, 1 < j ' < N — 1. 



We can confirm that 



e~ ux det 

l<j,k<N 



is the unique solution of this problem. 



The noncolliding Brownian motion with drift is defined as the system of Brownian motions 
conditioned that they never collide with each other forever. Then, the transition probability 
density is obtained by 

fit V\x) - hm ^(T-t,y) q ^t,y\x) 
p N {t,y\x) - T lim N u^ x) 

e-vy det [e"' Vk ] 

i<j,k<N v. { y-x)-MH/2 qN{tMx y 



e~ vx det \e v * Xk 

l<j,k<N L 



It is equal to (5). 
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Appendix B: A Combinatorial Limit of O'Connell Process 



Recently O'Connell introduced an interacting diffusive particle system [39]. Let 
i(>iP(x),x G M. N ,u G (the N- dimensional complex space) be the class-one Whittaker 
function, whose Givental integral representation is given by 



^\x) 



cxp 



(X) 
N-l j 



N 



£ Tj ' k ~ T 3-W 



dT, 



(Bl) 



_j=l \k=l k=l 

|e _ ( T j,*: _r 3+i,fc) _|_ e _ ( T j+i,fc+i _T j,fc)| 

3=1 k=l 

where the integral is performed over the space of all real lower triangular arrays with size 
N,T= (T jtk , 1 < k < j < N) conditioned T N>k = x k , 1 < k < N. The transition probability 
density P^ ,a of the O'Connell process with a parameter a > [40] is a unique solution of 
the equation 



d_ 
dt 



-A + W\og^\x/a) • V 



P%> a (t,y\x) = 



(B2) 



with the initial condition P^ ,a (0, y\x) = 5(x — y). The solution is given by 

Wv ( x / a ) 

with 



where i 



Q a N (t,y\x) = [ e-^l^hx/a^Ay/^s^dk, 
1 and sjv(a*) is the density function of the Sklyanin measure 



s N (n) = 



1 



(27T) 



[J |r(i(^-^-))l" 2 - 



(B3) 



(B4) 



(B5) 



l<j<^<Af 

We can show that (B2) is a geometric lifting with parameter a > of the diffusion 
equation (the backward Kolmogorov equation) of the noncolliding Brownian motion with 
drift [38, 40]. Then the BBO formula (5) is regarded as a combinatorial limit of the transition 
probability density of the O'Connell process in the sense that [51, 52] 



p%(t,y\x) = \imP£'' a (t,y\x), x, is G W,v, y G W N , t G [0,oo). 



(B6) 



We note that the above argument (the geometric lifting and combinatorial limit) gives 
another proof of the BBO formula (5) for the noncolliding Brownian motion with drift. 
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Appendix C: Geometric Brownian Motion 



For a^0,/i6l, and x > consider the linear stochastic differential equation (SDE) 

dX(t) = aX(t)dB(t) + iiX(t)dt, X(0) = x, t > 0, (Cf ) 

where B(t),t > is a standard Brownian motion started at 0. The parameters a and ji 
are called percentage volatility and percentage drift, respectively. The process X(t),t > 
is called a geometric (or exponential) Brownian motion [61]. The backward Kolmogorov 
equation for the process X(t) is given by 

did 2 d 

—u(t, x) = if 2 * 2 -^^! x ) + VXfaU(t, x). (C2) 

Set 

~ 1 ( ^ ^ 217+1 2 {r <l\ 

and put 

/ x \ a \ , x v ( ® 2 v 2 t (lny-lna;) 2 \ 
p(t;x,y) = -^=(xy) "exp — - 1 \^ } . (C4) 



2 V / 2^T ' V 2 2(7 2 * 

It is easy to confirm that (C4) satisfies the backward Kolmogorov equation (C2). Remark 

that p(t] x, y) is symmetric with respect to x and y; pit] x, y) = pit; y, x). 

Let 

m(dx) = —x^dx, (C5) 
a 2 

which is called the speed measure of X(t). By the general theory of one-dimensional diffusion 
processes [61], it is proved that for any A C K the probability P t (x,A) = Prob(X(t) G 
A|X(0) = x) is given by P t (x,A) = f A p(t;x,y)m(dy). In other words, if we set 

midy) 

P{t, v\x) = p{t; x, y)— — , (C6) 

dy 

then p(t, y\x) is the transition probability density of the process from x G R+ to y G En- 
during time t > 0; P t (x, -A) = J A p(t,y\x)dy. From (C4) and (C5), we have the following 



expressions, 



, . . 2 9 ~ ! IcH . . ~ / <j 2 z/ 2 t (lny — lnrr) 2 



<x 2 2 v / 2^t V 2 2 ^ 

1 /2/V / (? 2 v 2 t (my — lnrr) 2 



yy 

x) 



cxp 



y|a|V27rt Vx/ V 2 2a H 

yW hm e A-^ Hvlx) -~ vaH ) 2 }- (C7) 
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Let 

a t = \cr\Vi, 



Ht(x) = \nx + va 2 t = In a; + vo\. (C8) 

Then (C7) is written as 



1 f (my - fit(x)) 



p(t,y\x) = ^= ex PS Tn r- ( C9 ) 



It is nothing but the log- normal distribution with parameters a t and fit( x )- 
In the present paper we consider the case 



'2 



u = Q ^ ll =— (CIO) 



That is, the SDE is given by 



a 2 



dX(t) = aX(t)dB(t) + —X(t)dt, t > 0, (Cll) 

By a simple application of Ito's formula [70], we can confirm that the solution of this SDE 
is given by (25). 

Appendix D: Noncolliding Geometric Brownian Motion 

For arbitrary M G N, < t\ < • • • < t M < oo, the multitime joint probability density for 
the noncolliding Brownian motion without drift is given by [3, 14] 

p N (t 1 ,x W ;...;t M ,x (M) \x) 



M-l 



= n (4 M) -*r> ri ^w.-u^ L " N(t " x lT W T , odd 

l<j<k<N m=l lll<j<k<Ny x k ij) 

for any fixed initial configuration x = (x±, . . . , xn) G W^. 

By the transformation (25) from B(t) — > X(t),t > 0, the iV-particle system of geometric 
Brownian motions with percentage volatility a > conditioned never to collide with each 
other, which we call the noncolliding geometric Brownian motion, should have the multitime 
joint probability density in the form 

PN (tuy w ;---;tM,y^\y)= U (ln^-in^) 

l<j<k<N 

M-l gco/, \ 

II ggr(Wi-*m,y (m+1) ly (ro) ) n qN[ )\ V lV \ Y (D2) 
11 Yli< j<k <N( ln yk-^yj) 



X 

m=l 
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It does not seem to be possible to derive (21) from this general formula by just choosing an 
initial configuration y. 
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